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ABSTRACT 

In this Letter, the broad emission line (BEL) profiles of superluminal quasars 
with apparent jet velocities, (3 a > 10, (ultraluminal QSOs, or ULQSOs here- 
after) are studied as a diagnostic of the velocity field of the BEL emitting gas in 
quasars. The ULQSOs are useful because they satisfy a very strict kinematical 
constraint, their parsec scale jets must be propagating within 12° of the line of 
sight. We know the orientation of these objects with great certainty. The large 
BEL FWHM, ~ 3, OOOkm/s - 6, OOOkm/s, in ULQSOs tend to indicate that the 
BEL gas has a larger component of axial velocity (either random or in a wind) 
along the jet direction than previously thought. 

Subject headings: quasars: general — galaxies: active — accretion, accretion 
disks — black hole physics 



1. Introduction 



The most prominent feature of quasar optical/UV spectra are the conspicuous broad 
emission lines (BELs). The large line widths indicate a large spread in the velocity of the 
emitting gas. The nature of the velocity field of the BEL region is not well-constrained 
because the line of sight to the gas is not well defined in general. It is impossible to extricate 
the effects of orientation on the BEL FWHM (full width half maximum) from the effects of 
the central black hole mass and the luminosity of the accretion disk. The geometrical orien- 
tation of BLR (broad line region) and its relationship to the FWHM has received tremendous 
attention due to the potential consequences for studies of virial black hole estimates. 

Previous efforts to find orientation information on the BLR relied on one of the following 
methods: segregating fl at spectrum r adio loud quasars, FSQSRs, from steep spectrum radio 
lou d quasars, SSQSRs (ICorbin 111 9971) : correlating the FWHM with the radio spectral index, 



a (jJarvis and McClurd 120061 ); correlating FWHM with the core dominance parameter of 
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radio loud quasars, R (the re st frame flux density ratio at 5 GHz betwe en the unresolved core 
and re solved radio emission) . IWills and Browne I (119861 ); or with Ry of IWills and Brotherton 
( 119951 ) (the ratio of the radio core flux density at 5 GHz to the optical flux density which is 
isotropic, if it is from an accretion disk); the rotation of the polarization vector across the 
BEL can constrain the scattering geometry has been applied to some Seyfert 1 galaxies and 
supports a disk-like flattened geometry of the BLR (ISmith et alll2005l ). 

These methods give vague orientation information which has limited our ability to re- 
solve the geometry of the BEL velocity field. In particular, SSQSRs are probably concen- 
trated in the range of 9 < 55° lines of sigh t from the jet axis based on the statistics of 
the complete 3CRR sample (jOgle et all 120061 ). Due to relativistic Doppler boosting, the flat 
spectrum radio core flux, if viewed at 9 < 20°, will be drastically enha nced so as to swamp 
the s teep lobe flux at 5 GHz and the source will appear as a FSQSR (ILind and Blandford 
1985h . Thus, for SSQSRS, we expect 20° < 9 < 55°. The FSQSRs, are restricted to 
lines of sight within 9 = 30° from the jet axis. This constraint on the FSQSRs is much 
weaker than is generally acknowledged and provides tremendous overlap between the two 
sub-populations. We illustrate this with an example. Recall that for an unresolved source 
the observed flux density, S u (obs), is Doppler enhanced relative to the intrinsic luminosity, 
S u (obs) = 5 3+a S l ,(intr), where the Doppler factor, S is given in terms of T, the Lorentz factor 
of the outflow; cj3 = sjv* + v y + v'j, the three velocity o f the outflow; and a, is the spectral 
index: 5 = 1/[T(1 — /3cos9)} (ILind and Blandfordl Il985l ). Consider an FR II quasar that is 
viewed at 9 = 50°. It is a "typical" powerful FR II quasar at z=0.5, with steep spectrum 
lobes S u (obs)=3 Jy at 178 MHz and the lobe has a — 1. The radio core is chosen to have 25 
mJy at 5 GHz with a = 0. Then a = 0.85 between 2.7 GHz and 5 GHz, a very typical value 
for a SSQSR. If the jet is moving with (3 = 0.995, (r = 10) and the same object were viewed 
at 9 = 30° then the observed core flux density is about 500 mJy at 5 GHz and the lobe 
flux is only about 100 mJy at 5 GHz. The same object that was a typical lobe dominated 
SSQSR viewed at 9 = 50° is transformed into a typical FSQSR (a = 0.23 between 2.7 GHz 
and 5 GHz) if it is viewed at 9 = 30°. This example indicates that FSQSRs should exist 
for lines of sight 9 = 30° from the jet axis which overlaps the range of 9 for SSQSRs (based 
on the angular range discussed above, the average value for the lines of sight to SSQSRs is 
9 = 39°)The difference between 30° and 12° (which is the maximum 9 of the ULQSO sample 
defined below) is highly significant for this analysis because the projection of an equatorial 
velocity component onto the line of sight ~ sin^ (c.f. equations (3.1) and (3.2)): this equates 
to a factor of 2.5 for the two values of 9. It is concluded that segregating FR II QSOs by 
spectral index gives very crude orientation information. 



The quantities R and Ry give improved angular information because it captures the 
Doppler enhancement factor. However, this is offset by the fact that intrinsic values of 
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radio core luminosity, radio lobe luminosity and accretion disk luminosity are not "stan- 
dard candles," but vary considerably from quasar to quasar. Secondly, there is large opti- 
cal contamination of the accretion disk luminosity from the Doppler boosted jet in many 
FSQ SRs that further s mears t he angular resolution o f this method. The correlation analysis 



of IWills and Browne I (119861 ); IWills and Brothertonl (119951 ) is convincing, but it really only 



tells us that there is some component of the BLR that is concentrated near the equatorial 
plane, it doesn't tell us if this represents the bulk of the BEL emitting gas or a much smaller 
fraction like 1/3 of the total BEL gas (cf. eqn (2.2)), the correlation could exist in either 
case. 

Our understanding of the velocity field of the BEL gas has been limited by the poor 
orientation resolution of these methods. In the end, one is left with a lar ge unc e rtaint y 



in the orientation factors that are relevant for virial black hole estimates iKrolikl (120011 ). 
Thusly motivated, in this paper, we study a sample of objects in which the line of sight is 
very well-constrained. Even though the sample is not large, 76 broad lines from 33 QSOs, 
the accurate orientation information is very insightful. The sample consists of ultraluminal 
quasars (ULQSOs, hereafter) that are defined by the condition that their projected apparent 
velocity on the plane of the sky is, (3 a > 10. Since, f3 a = (/? sin#)/(l — j3 cos 9), the (3 a > 10 
condition and the fact the intrinsic velocity of the jet plasma is subluminal (3 < 1, means 
that the line of sight makes an angle 9 < 12° to the jet axis. With the knowledge that we are 
viewing the quasars extremely close to the polar axis, implies that the FWHM of the BELs 
will give unambiguous information on Doppler shift associated with one component of the 
velocity field of the BEL emitting gas, the axial component along the jet axis. This article 
considers all of the prominent BELs that were accessible to observation, H/3, Mgll, CIII] and 
CIV. There is strong blending between CIII] and Silll, thus the analysis of these lines is not 
treated with as much weight as the other BELs in the following. The FWHM that are found, 
correspond to velocities that average m 4,000 km/sec for all BELs. The implication is that 
QSOs have prodigious quantities of high velocity gas moving both parallel and antiparallel 
to lines of sight along the jet direction. Formally, the results of this paper are restricted to 
radio loud AGN. The extent to which it applies to radio quiet AGN is unclear. 



2. The Ultraluminal Quasar and Comparison Samples 

Fortunately, there are three large multi-epoch VLBI studies. The sample in table 1 



Jorstad et al 



repres e nts all the (3„. > 10 AGN that were found by iKellermann et al I (120041 ) ; 
(120011 ); lYermeulen and Cohen! (119941 ) for which there is some archival FWH M data. There 



are many preferences for de-constructing BELs into various components (IMarziani et al 
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20031 ; IVestergaard and Peterson! 120061 ) . There is likely much merit in this endeavor, but the 
decompositions and FWHM estimates are not produced uniformly in the references cited 
in tablel and table 2 (described below). The choice is to either greatly limit the sample 
size to a few sources from one reference or introduce scatter depending on whether narrow 
components were subtracted or how the very broad component was fit. A very small sample 
of ULQSOs is not of much interest from a statistical point of view. The FWHM in the 
control samples in table 2 are drawn from largely the same references that appear in table 
1. Thus, the degree of scatter from different derivation techniques occurs in all samples and 
is likely to be quantitatively similar. In order to minimize this scatter, an effort was made 
to find single Gaussian fits to the data, if they were available. This gives us the maximum 
amount of commonality with the control samples in table 2 f or the low ionization B ELs 
because man y of these FWHM com e from single Ga ussian fits: lLawrence et al.l (1996). for 
the ULQSOs; lLawrence et al.l dl99d): ICorbin I Jl997h for the sample of iGu et all (bpOlh : and 



Aars et all (boosh 



Lawrence et al 



<jl996h for the 3CRR sample. It would be important to 
followup this work with a uniform reanalysis of the raw data that segregates the various 
components for comparison. This is far beyond the scope of this effort. 

In order to obtain more information than just the fact that there are large BEL velocities 
along the jet axis, it is of interest to compare this data with samples that represent sources 
viewed from far more oblique angles. Unfortunately, we can not use the large database of ra- 
dio quiet quasars because there are numerous claims in the literature that the BELs of radio 
quiet and radio loud quasars have significant d ifferences. For example, som e authors claim 
that B./3 FWHM is larger in radio loud QSOs (IMcClure and Dunlopll200ll ). Thus, we must 
restrict our discussion to radio loud objects and the SSRQRs are the sources that are likely 
viewed at relatively steep angles to the line of sight. Given that we are dealing with a rela- 
tively small number of sources ~ 20 for a particular ULQSO BEL and samples of SSQSRs in 
the literature are of modest size ~ 50, it is unrealistic to match all parameters (i.e., redshift, 
UV luminosity, extended radio power). First of all, it is impossible to match UV power ap- 
propriately because the accretion disk emission is contaminated by the jet optical emission in 
the ULQSOs. It is also extremely difficult to match extended radio power, since ULQSOs are 
core dominated and one needs very deep, high dynamic range radio maps (to find the rela- 
tively faint extended flux) which do not exist except for a few of the ULQSOs. Thus, we chose 
to match the samples by redshift since this is often related to one or more of these other pa- 
rameters by selection biases. The samples used for comparison and contrast are listed in table 
2 and the distributions of FWHM can be directly compared. As a check of our sample selec - 
tion method, we can compare to the carefully selected samples of iJarvis and McClurd (120061 ). 
They find for a combination of H/3 and Mgll, FWHM(FSQSR)/FWHM(SSQSR)=0.77, and 
we find for H/3; FWHM(ULQSO)/FWHM(SSQSR)=0.67, for the combined SSQSR sample 
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of iGu et all (1200 ll ); ICorbin I (Il997l ). which is consistent with the ULQSOs being viewed a 
little closer to the pole than a typical FSQSR as discussed in the introduction and verifies 
the lack of a significant selection bias to our study. 



First we compa re H/? to the sample of SSQSRS from IGu et all (120011 ) as was done in 
Oshlack et all (120021 ) in their comparison with FSQSRS. The sample is comprised of all the 
SSQSRS with > 1 Jy at 5 GHz from the S4 and S5 surveys and is 97% complete (IGu et al 
200ll ). This is well matched in redshift to the ULQSOs, but it is small. So, we construct 
a larger H/? control sample by adding the ICorbin I (119971 ) sample of SSQSRS which is at 
similar redshift. The essential conclusions of this article are unaffected by which sample 
is used as the control. There was no existing list of Mgll FWHM for SSQSRs, so for a 
comparison the noncompact (> 20 kpc ) SSQSRs from the 3CRR catalog that are known 
not to be dominated by a powerful blazar core were chosen to represent the steep lines of 
sight. T he samples for CHI] and CIV are drawn from the same references, the high z sources 
ar e fromlCorbin and Francis I (11994 ) and the low redshift contribution is from the HST data 
of IWills et all Jl995h in thetables. 



From table 2, based on the size of the standard deviation, from a statistical point of view 
one can not claim that there are differences between the SSQSR FWHMs and the ULQSO 
FWHMs for any of the four BELs. However, we note that the means and medians for H/? 
and Mgll FWHMs are smaller for the ULQSOs and the opposite is true for CIII] and CIV. 
We consider the popular notion that the BEL gas is distributed in an equatorial pancake, 
orthogonal to the jet axis, with a random velocity, v r , superimpo sed on an equatorial v elocit y, 



v p , that is predominantly bulk motion from Keplerian rotation iJarvis and McClurd (120061 ). 



FWHM 



+ vl sin 2 9 . 



In lLabita et all (120061 ). they assume that the random component is from Keplerian orbits of 
random inclination and find the relation 



FWHM = 2(ci-^= sin0 + c 2 -j=)v B lr 



(2-2) 



where C\ and c 2 "represent the relative fraction of the FWHM due to the" planar and isotropic 
components, respectively. Using the angular coordinate ranges motivated in the Introduction 
(0° < 9 < 12° for ULQSOs and 20° < 9 < 55° for SSQSRs) and the mean FWHM from 
Table 2, one can estimate the ratio v r /v p for each BEL after appropriately averaging (3.1). 
For example, for H/?, we solve the equation, 



12 c 



cos 0° - cos 12° 

/3898\ / 1 
" V5836/ Vcos20°-cos55° 



v 2 + sin 9 sin 9 d9 



55° 



20-- 



v 2 + v 2 sin 9 sin 9 d9 . 



(2-3) 
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This equation can be solved numerically to give, v r /v p = 0.537 for H/3. Similarly, for Mg II, 
v r /v p = 0.502 and for CIV (and CIII), the only consistent solution is v p = 0. If the mean 
values of the FWHM from table 2 are replaced in (3.3) by the median values, the results 
are slightly changed; for H/3, v r /v p = 0.563; for Mg II, v r /v p = 0.496; and for CIV (and 
CIII), the only consistent solution is v p = 0. This makes the result seem more robust. If 
one assumes that v r and v p represent the rms velocity of the two BE L components then one 
can cl aim approximate equality in (2.1), as in iMcClure and Dunlopl (1200 ll ); lLamastra et al 
( 120061 ). Using the mean values in table 2 transforms the average of (2.1) over solid angle in 
(2.3) into two equations in two unknowns and one can solve for the mean values of v r and 
v p : v r = 1877 km/ s, v p = 3496 km/s for H/3; v r = 1979 km/ s, v p = 3942 km/s for Mgll; 
and v r = 2150 km/ s, v p ~ km/s for CIV. Similarly, one can average (2.2) over solid angle, 
to find cx/c 2 = 1.04 for H/3 and Ci/c 2 = 0.93 for Mgll and c x jc 2 = for CIV. 

The improved angular resolution of this study can be compared to previous analysis 
of the BLR geometry. From (2.1) and (2.2), the low ionization BLR can be characterized 
by a mixture of planar and random components. For CIV, the model breaks down and 
the only consistent solution is purely random motion within the model. The value of v r 
is large r than has what ha s been estimated by lower angular resolution studies. I n section 



II b o f IWills and Browne I (119861 ). it is estimated that 0.15 < v r /v p < 0.31 and in ICorbin 



( 119971 ) v r /vp ~ 0.15 for H/3, compared to our value of 0.537. Similarly, iNikolajuk et all (120061 ) 
assume that (v r /v p ) 2 <C 0.28, while we find to the contrary that (v r /v p ) 2 ~ 0.3. Furthermore 



the ULQSO analysis applied to CIV strongly contradicts the findings of lLabita et all (120061 ) 
who assert that the CIV FWHM in quasars are consistent with a planar distribution of BEL 
gas near the equator. The discrepancy could be a manifestation of the difference between 
their radio quiet sample and our powerful jet sample. This effort seems to be the first that 
actually performs the numerical averaging of (2.1). 



3. Conclusion 

In this paper, we study the FWHM of the BELs of ULQSOs. We applied our findings 
to the planar models of the BLR, focusing on the velocity components of the BEL gas 
both redward and blueward along the jet axis (perpendicular to the inner accretion disk). 
Our analysis of the ULQSOs supports the notion that the high ionization BELs come from 
gas with a different velocity field than the low ionization BELs. The ULQSOs FWHM 
tends to indicate that the putative planar flow has a much larger random component than 
previously acknowledged and the planarity of the low ionization BLR has been overstated 
in the literature, since it was shown in section 3 that C\jci ~ 1 in equation (2.2) (from the 
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methods of [Labita et all (120061 ) this implies that the isotropic gas and the planar flow are 
equal contributors to the FWHM). Formally, these results apply only to radio loud quasars. 
One might be suspicious that the results presented are skewed by sample selection biases. 
First of all, the anomalously large axial velocity component in the context of the equatorial 
low ionization BLR model is corroborated by three control samples, 2 of which are virtually 
complete. There does not appear to be any obvious selection bias in choosing ULQSOs. It 
was noted in section 3 that a selection bi as does not appear to exist b ased on the consistency 



with the carefully matched samples in iJarvis and McClurd (120061 ). The large redshifted 



BEL velocities in ULQSOs is concluded to be robust and the dynamics of this phenomenon 
requires further theoretical investigation. Any deductions beyond this are speculative due 
to the large scatter in the data. 
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Table 1: Broad Line Widths of intraluminal Quasars 



Source 


z 


max /3 a 


CIV FWHM 


CIII] FWHM 


Mgll FWHM 


H/3 FWHM 


ref b 








(km/s) 


(km/s) 


(km/s) 


(km/s) 




0016+731 


1.781 


11.9 


5610 


5390 


4172 




1, 4 


0106+013 


2.107 


23.3 


5100 


5700 






1, 5 


0332+321 


1.263 


24.5 


4000 


5700 


2900 




1, 6 


0336-019 


0.852 


12.4 








4875 


2, 7 


0420-014 


0.915 


14.2 








3480 


1, 8 


0736+017 


0.191 


11.5 








4628 


1, 7 


0827+243 


0.939 


27 


6099 








2, 9 


0836+710 


2.16 


14.1 


6775 


10403 


3047 


3410 


2, 4, 10 


0850+581 


1.322 


12.7 




5107 


6407 




1, 4 


0859-140 


1.327 


16.3 


4500 


5400 


4400 


5700 


1, 6, 11 


0906+015 


1.108 


11.8 




6160 


3990 




1, 8 


0945+408 


1.252 


22.5 


3454 


6003 


7593 




1, 4 


0953+254 


0.712 


12.4 








3990 


1, 7 


1055+201 


1.11 


10 


3698 


2979 






1, 9 


1127-145 


1.187 


19.8 


4880 


6190 


2750 




2, 8 


1156+295 


0.729 


14.1 


5072 


5601 


2664 


3700 


1, 9, 12 


1226+023 


0.158 


14.1 


3941 






4040 


2, 13, 14 


1253-055 


0.536 


10.4 


5502 


2769 


4037 


3100 


3, 15, 16 


1308+326 


0.996 


16.3 






4016 




1, 17 


1406-076 


1.493 


28.4 




5800 


3650 




2, 8 


1508-055 


1.191 


31.4 


1680 


2160 


2230 




1, 8 


1510-089 


0.36 


18.9 






2680 


2796 


1, 8, 13 


1611+343 


1.401 


24 


4744 


5575 




5600 


1, 11, 16 


1633+382 


1.814 


11.5 


4332 


4840 


9650 




2, 4 


1641+399 


0.594 


17 


4829 


4098 


4196 


4320 


1, 4, 16 


1642+690 


0.751 


16 




3190 


4896 


1850 


1,4 


1655+077 


0.621 


15.3 






2080 




1, 8 


1828+487 


0.691 


14.8 






3010 


4245 


1, 4, 18 


1901+319 


0.635 


15.2 








2580 


2, 19 


2223-052 


1.404 


32.5 


3449 


4418 






1, 16 


2243-123 


0.63 


10.7 








4258 


1, 18 


2230+114 


1.032 


11.7 


3340 


3903 






2, 16 


2251+158 


0.859 


15.7 


3145 


5175 




3700 


2, 7, 16 



"Column (1) is the source, followed by the redshift and the maximum f3 a in the three studies. Columns 
(4), (5), (6) and (7) are the FWHM of CIV, CIII], Mgll and H/3, respectively. Column (8) is a list of the 
references, the first reference is to the VLBI data and the remainder are for the FWHM. 



b Rcfcrences: 


1. Kellermann et al (2004) 2. Jorstad et al (2001) 3. Vermeulen and Cohen (1994) 4. 


Lawrence et al. (1996) 5. 


Corbin and Francis (11994) 6. Brotherton et al (19941) 7. Jackson and Browne! 


(1991) 8. Wills and Browne (1986 


) 9. Kuraszkiewicz et al (2004) 10. Mcintosh et al (1999) 11. 


Netzer et all ( 


1995 s ) 12. Brothertonl ( 


1996) 13. Oshlack et all (2002) 14. Vestereaard and Peterson (2006) 


15. Wilkes 


dl986h 16. 


Wills et al. (1995) 17. Wans et al. (|2004) 18. Marziani et al ()2003h 19. 


Gelderman and Whittle ( '. 


99l 
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Table 2: Broad Line Widths of intraluminal Quasars versus Steep Spectrum Radio Loud QSOs a 



Line 




ULQSOs 










SSQSRs 








FWHM 


median 


z 


number 


FWHM 


median 


z 


number 


ref b 




(km/s) 








(km/s) 










H/3 


3898 


3990 


0.80 


17 


5460 


4500 


0.59 


29 


1 




±1011 




±0.50 




±3639 




±0.39 






H/3 


3898 


3990 


0.80 


17 


5836 


5530 


0.51 


52 


1, 2 




±1011 




±0.50 




±3476 




±0.32 






Mgll 


4125 


3990 


1.11 


19 


6398 


6222 


0.98 


24 


3CRR 3, 4, 5, 6 




±1997 




±0.52 




±2668 




±0.39 






CIII] 


5074 


5390 


1.26 


21 


3904 


3700 


1.59 


43 


6, 7 




±1707 




±0.45 




±1698 




±0.68 






CIV 


4429 


4500 


1.27 


19 


4279 


3900 


1.60 


45 


6, 7 




±1189 




±0.47 




±1369 




±0.68 







a The first column in table 2 is the BEL under consideration. Next, we list the ULQSO mean and median 
values of the FWHMs from table 1 followed by the number of objects in the subsample in columns (2)-(5). 
Columns (6)-(9) are the same for the samples of SSQSRs. 

( 1996h 4. lAars et all (|2005h 5. 



b References: 1. 


Gu et al 


Brotherton et al 


1994) 6. 



Corbin (1997) 3. 



Lawrence et al 



